Statistical mechanics of Semiflexible Bundles of Wormlike Chains 
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We demonstrate that a semiflexible bundle of wormlike chains exhibits a state-dependent bending 
stiffness that alters fundamentally its scaling behavior with respect to the standard wormlike chain. 
We explore the equilibrium conformational and mechanical behavior of wormlike bundles in isolation, 
in crosslinked networks, and in solution. 
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In recent decades, the wormlike chain (WLC) has 
emerged as the standard model for the description of 
semiflexible polymers [1]. The defining property of a 
WLC is a mechanical bending stiffness, k/, that is 
an intrinsic material constant of the polymer. Within 
this framework, numerous correlation and response func- 
tions have been calculated, providing a comprehensive 
picture of the equilibrium and dynamical properties of 
WLCs [3,i,ll. A number of experimental studies have 
demonstrated the applicability of the WLC model to 
DNA 3 and F-actin [6], among other biological and 
synthetic polymers. Significant progress has also been 
made towards the description of the collective properties 
of WLCs, for example, in the form of entangled solu- 
tions. One of the hallmarks of this development is the 
scaling of theplateau shear modulus with concentration, 
which is well established experimen- 
tally 

Another important emerging class of semiflexible poly- 
mers consists of bundles of WLCs 
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Semiflexible 

polymer bundles consisting of F-actin or microtubules 
are ubiquitous in biology 1^, and have unique mechan- 
ical properties that may well be exploited in the design 

IS [3 



of nanomaterials [13| . As shown by Bathe et 
wormlike bundles (WLB) have a state-dependent bend- 
ing stiffness, kb, that derives from a generic interplay 
between the high stiffness of individual fllaments and 
their rather soft relative sliding motion. In this Letter, 
we demonstrate that this state-dependence gives rise to 
fundamentally new behavior that cannot be reproduced 
trivially using existing relations for WLCs. We explore 
the consequences of a state-dependent bending stiffness 
on the statistical mechanics of isolated WLBs, as well as 
on the scaling behavior of their entangled solutions and 
crosslinked networks. 

We consider the bending of ordered bundles with 
isotropic cross-section. A bundle consists of N filaments 
of length L and bending stiffness Kf. Filaments are ir- 
reversibly crosslinked to their nearest neighbors by dis- 
crete crosslinks with mean axial spacing 5. Crosslinks 
are modeled to be compliant in shear along the bundle 
axis with finite shear stiffness /cx, and to be inextensible 



transverse to the bundle axis, thus fixing the interfila- 
ment distance, b [i^. Bundle deformations are charac- 
terized by the transverse deflection r^(s) of the bundle 
neutral surface at axial position s along the backbone 
and by the stretching deformation Ui{s) of filament i. 
The torsional stiffness of the bundle is assumed to be of 
the same order as the bending stiffness. Thus, as long as 
transverse deflections remain small ("weakly bending") 
the two components of rj^ are decoupled and effects of 
twist are of higher order 19]. The bundle response may 
then be analyzed in planar deformation, where the bend- 
ing stiffness results from the superposition of 2M = ^/N 
bundle layers. 

The WLB Hamiltonian consists of three contributions, 
Hw-LB = Hhcnd + -^stretch + -^shear- The first term corre- 
sponds to the standard WLC Hamiltonian 
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which is the same for each of the iV filaments. The second 
term accounts for filament stretching. 



stretch 



Mk.S / ds 
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where kg is the single filament stretching stiffness on 
the scale of the crosslink spacing 5. No particular form 
for bending and stretching stiffnesses is assumed, but 
one may think of the filaments as homogeneous elastic 
beams with Youngs modulus i?, for which Kf Eb'^ and 
ks ^ Eb'^ /S. Alternatively, ks may represent the entropic 
elasticity of a WLC, for which kg ^ n'j/TS'^. 

The third energy contribution, iJshcar, results from the 
crosslink-induced coupling of neighboring filaments. To 
minimize the crosslink energy, any relative filament slip 
induced by cross-sectional rotations 9 — dsr± = r'j_ must 
be compensated by filament stretching (Fig. [1]). This 
crosslink shear energy, which simply suppresses relative 
sliding motion of neighboring filaments, is given by 
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FIG. 1: Illustration of the geometry of a single bundle layer 
(The full bundle consists of 2M layers that are stacked in 
parallel.). The bundle is deflected through the angle 6 = r'^. 
If filament i stretches the amount m = Ui+i+b9, the crosslink 
(dashed line) remains undeformed with zero shear energy. 



where Aui — Ui — Ui-i. A related model for two filaments 
was introduced by Everaers et al. 



in Ref. [17|, where 
special emphasis was placed on the limit of inextensible 
filaments, kg — > oo. In that model, the anisotropic bundle 
cross-section leads to a coupling of in-plane and out-of- 
plane bending modes [3| that is absent in the present 
model because it has a symmetric cross-section. 

Functional differentiation of the Hamiltonian results in 
the (overdamped) equations of motion 



N^frT-^^Y.(^< + brl) = F(r^,.), (4) 

i 

ksSu'^ + ^iAu,+i- Au,) = 0, (5) 





where is a transverse force that may represent fluid 
drag, random thermal noise, or other external loading. 
To proceed, Eq. ([5]) is solved together with appropri- 
ate boundary conditions, so as to eliminate the Ui in 
Eq. (|4]). The calculations are most easily performed in 
Fourier-space, where we write for the expansions r±{s) — 
J2n'^nSm{mTs/L) and Ui{s) = X]„ "^n cos(n7rs/L), ap- 
plicable to pinned boundary conditions. The resulting 
equation of motion for r„ then takes the simple form 
i^nQni^n = Frn with a modc- number dependent effective 
bending stiffness k„. The general result for k„ is ob- 
tained using the standard ansatz ~ w*, which reduces 
Eq. ([5]) to an equation that is quadratic in w. 

In the following, we present an approximate solution 
to Eqs. (Ill) and ^ that is based on the assumption that 
filament stretching increases linearly through the bundle 
cross-section, Ui ~ Au • (i + 1/2) 27]. Although compar- 
ison with the exact solution demonstrates that Ui in gen- 
eral varies nonlinearly with i [20.] . it turns out that the 
effective bending stiffness k„ is insensitive to this non- 
linearity. At the same time, the linearization simplifies 
the formulas substantially, so that the effective bending 



stiffness is given in closed form by 
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Kn = Nkj 
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with a dimensionless bending stiffness k 



f 



(6) 



and a length-scale A = [L/ ^)^J Mkf / {M - 1/2), that 
depends on the shear stiffness ky, via the dimensionless 
coupling parameter a = kxL'^ /ksS'^ . 

For any given mode- number ~ n/L, three differ- 
ent elastic regimes emerge as asymptotic solutions for 

» 1 and respective values of a [l^[l^. For large shear 
stiffness {a ^ N), the fully coupled bending scenario is 
obtained, where the bundle behaves like a homogeneous 
beam with k„ ^ N'^kg. For intermediate values of the 
shear stiffness (1 ^ a ^ N), the bending stiffness in 
the shear dominated regime is k„ ~ NkxQn'^ and the full 
mode-number dependence of Eq. ^ has to be accounted 
for. Finally, decoupled bending of N laterally indepen- 
dent, but transversly constrained, filaments is found in 
the limit of small cross- link shear stiffness (a <C 1), where 
the bending stiffness is simply k„ = Nnf. 

In the particular limit of N —^ oo and fixed bundle 
diameter D — b\/N ^ L, Eq. ( Et r educes to the Timo- 
shenko model for beam bending [21], which was recently 
used to interpret bending stiffness measurements on mi- 
crotubules 22, 2^ and carbon nanotube bundles [l^. In 
this limit 
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where we have used the expressions of kg and k f for ho- 
mogeneous beams and defined G = ky,/d. While this 
limit serves as a consistency check for our mathemati- 
cal analysis, real bundles consist of a finite, and often 
small, number of constituent filaments, for which Eq. ^ 
cannot be applied to describe the full range of bending 
behavior captured by Eq. ([6]). Indeed, in Eq. ([7]) no 
decoupled bending regime exists and the bending stiff- 
ness vanishes as the crosslink shear stiffness approaches 
zero |29|]. The condition, a ^ N, delineating the remain- 
ing two regimes can be rewritten as E/G ~ {L/D)'^ ^ 1, 
which re-emphasizes the small value of cross-link shear 
stiffness in the intermediate regime. 

For fixed values of {N,a), the bundle bending stiff- 
ness Eq. ^ crosses over from fully coupled to decoupled 
bending via the intermediate regime as the mode-number 
Qn is increased. Thus, different modes may belong to dif- 
ferent elastic regimes, rendering the fluctuation proper- 
ties of the bundle non-trivial and qualitatively different 
from single semiflexible polymers. This cross-over is me- 
diated by the length-scale A, which acts as a cut-off on 
the fluctuation spectrum: whereas wavelengths 'ti X 
belonging to the decoupled regime are characterized by a 
constant bending stiffness, modes with 3> A acquire 
a higher stiffness k„ ~ q~'^ and are thereby suppressed. 
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Finally, for even longer wavelengths ^ A\/lV, the 
bending stiffness reattains a constant, limiting value. As 
an example (taken from Ref. [l3|) we found A ~ for 
actin/fascin bundles with N « 30, L « 50/iTO. 

In situations where modes pertaining to the interme- 
diate regime are irrelevant, the g-dependence of k„ drops 
out and one recovers the single WLC result, albeit with 
a renormalized persistence length Ip Nip in the de- 
coupled, and Ip — > N'^lp in the fully coupled, regimes, 
respectively. In other cases, calculation of the tangent- 
tangent correlation function demonstrates that the per- 
sistence length cannot be defined unambiguously. As in- 
dicated in Ref. the correlation function does not 
decay exponentially, but rather exhibits a complex struc- 
ture at intermediate distances [23|. In the following, we 
will therefore explore the consequences on the statisti- 
cal mechnanics of the WLB in particular as regards the 
intermediate regime. 

First, consider the force-extension relation as cal- 
culated from the end-to-end distance R{F) = L — 
f^sT/ {Kriai + F), where F is the force applied to the 
bundle ends [30|. For small stretching forces one may 
readily calculate the linear response coefficient fcontr — 
F/{R{F) — i?(0)) using a Taylor series expansion. The 
result in the intermediate regime is 

« r^3?C ' » ^ » ^) • (8) 

which is inversely proportional to bundle length, like a 
mechanical beam. Importantly, the strong dependence 
of fcentr(i) ~ applicable to single filaments (and 

the other two regimes) is lost. This has dramatic con- 
sequences on the plateau value of the shear modulus 
G in crosslinked bundle networks, which in affine the- 
ories [24] is assumed to be given in terms of /contr by 
G ^ fcontr(0/'?j where the mesh-size ^ depends on con- 
centration c as ^ ~ c~^^'^. Accordingly, in the interme- 
diate regime one finds G ^ c, which is a much weaker 
concentration-dependence than G ~ c^/^ [i^ applica- 
ble to single filaments. It is worthwhile noting that the 
force-extension relation is strongly nonlinear (see Figl2|) , 
rendering the linear response valid only for very small 
relative extensions. In this particular example the lin- 
ear response formula deviates from the exact solution by 
50% at only « 3% and « 0.7% strain in the decoupled 
and the fully coupled limits, respectively. 

Bundle behaviour under compressive forces further 
highlights the unusual properties of WLBs. Because the 
bending stiffness in the intermediate regime scales with 
the length of the bundle a,s kb ^ L^, the Euler buckling 
force Fc ~ kb/L"^ ^ Nnf/X^ is independent of bundle 
length. This unique property may well be exploited in 
polymerizing biological bundles such as filopodia, which 
may increase their contour length against compressive 
loads without loss of mechanical stability. 




FIG. 2: End-to-end distance R{F)/L as a function of stretch- 
ing force FL^ / Kf for a bundle oi N = A filaments and L = Ip. 
The black curves correspond to X/L — 0.01, 0.1, ..,0.7. Thick 
red curves relate to (bottom) decoupled and (top) fully cou- 
pled bending, respectively. Dashed lines correspond to the 
respective linear response regimes. 

Complementary to the elasticity of crosslinked net- 
works of WLBs, we turn next to the elasticity of their 
entangled solutions. The generally accepted theory for 
the concentration dependence of the plateau modulus of 
entangled WLCs is based on the free energy change AF 
of confining a polymer to a tube of diameter d [3, . 
The associated change in free energy is written as AF ~ 
ksTL/ld, which defines the deflection length Id to be the 
scale at which the polymer starts to interact with its en- 
closing tube. The deflection length itself is connected to 
the tube diameter d and the filament concentration c via 
the standard excluded volume argument fo*], Qd = ki/cL, 
which balances the excluded volume of the tube with the 
available volume per filament. All that remains is the 
calculation of the tube diameter c? of a single polymer 
confined by the potential 

where the confinement length Ic is defined as a measure 
of the strength of the potential. While Ic = U in the 
standard WLC, we will see shortly that this does not hold 
for WLBs in the intermediate regime. First, consider 
the transverse fluctuations of an unconfined bundle, in 
particular the average value d^ = /s(''J-(5)^)- This is 
most easily calculated as 

dl - LX^/Nlp , {XVn > L > A) (10) 

which has to be compared to the WLC result for which 
dp ~ /Ip- In the presence of the confining potential, 
the same calculation yields 

d^^ljX/Nlp, (AViV>Zc>A). (11) 
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FIG. 3: Tube diameter /{l^/Nlp) as a function of con- 
tour length L/lc for various X/lc and M = 20. Thick (red) 
curves correspond to (top) decoupled and (bottom) fully cou- 
pled bending, respectively. For short filaments the interme- 
diate regime is visible through the linear slope d? ^ L (see 
Eq. (Unj). For long filaments the fluctuations saturate. By 
increasing A the tube is becoming wider (Eq. 



For strong confinement Ic ^ A, the potential suppresses 
all modes of the intermediate regime and one recovers the 
expression valid for single filaments, <P ~ l^/lp- The gen- 
eral result for the tube diameter is depicted in Fig. [31 As 
the contour length L of the bundle is increased it begins 
to "feel" the presence of its enclosing tube at the deflec- 
tion length L = Id- By comparing Eq. PU)) with Eq. (fTTj) 
one finds Id ~ ^c/-^' which is valid in the intermediate 
regime. At the same time, Id = Ic in the decoupled and 
fully coupled regimes, where the deflection and confine- 
ment lengths are identical. 

One may use these results to rewrite the deflection 
length as a function of concentration c. In the intermedi- 
ate regime the result is 1^ ^ Nip/ [XcL)'^ , which replaces 
the usual result ~ Nip/ {cL}^ valid in the decoupled 
regime (strong confinement). The free energy of confine- 
ment and the elastic plateau modulus G ^ {cL)AF/L 
now depend on A and thus on the properties and density 
of the crosslinks. The modulus displays a cross-over that 
is mediated by concentration, 



G - ksT 



'{cL)^/HNlp)-^/^X^/\ c<^c\ 
{cLy/^{Nlp)-^/\ c^c*, 



(12) 



where we defined the cross-over concentration as (cL)* ^ 
y/]v7^A"5/2. Below the even smaller concentration c** ^ 
c*N-^/^, the fully coupled regime is entered and the 
modulus again scales as G ~ c"^^^ . 

Having addressed equilibrium properties of WLBs, fur- 
ther consequences of the state-dependent bending stiff- 
ness on dynamic response functions remain to be ex- 
plored, along with the effects of nonperrnanent crosslinks. 
Additional experiments 



and synthetic bundles. 
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